Self-assembly of multicomponent structures in and out of equilibrium 
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Theories of phase change and self-assembly often invoke the idea of a 'quasiequilibrium', a regime 
in which the nonequilibrium association of building blocks results nonetheless in a structure whose 
properties are determined solely by an underlying free energy landscape. Here we study a proto- 
typical example of multicomponent self-assembly, a one-dimensional fiber grown from red and blue 
blocks. If the equilibrium structure possesses compositional correlations different from those charac- 
teristic of random mixing, then it cannot be generated without error at any finite growth rate: there 
is no quasiequilibrium regime. However, by exploiting dynamic scaling, structures characteristic of 
equilibrium at one point in phase space can be generated, without error, arbitrarily far from equi- 
librium. Our results thus suggest a 'nonperturbative' strategy for multicomponent self-assembly in 
which the target structure is, by design, not the equilibrium one. 



Many theories of phase change and self-assembly place 
at their heart the idea that dynamical trajectories follow 
low-lying paths on the free energy landscape connect- 
ing reactants and products p~HTTj . This idea underpins 
rate equation theories [1, 3-5 , classical nucleation the- 
ory [TJ [J and density functional theory [7HTT] (which 
assume that a structure's morphology is determined by 
minima or low-lying paths on the underlying free en- 
ergy landscape), and the conjecture of Stranski and To- 
tomanow [12], which states that a system, confronted 
by a set of free energy barriers, will evolve by crossing 
the lowest of them. These formal statements reflect the 
intuition that one can generate structures characteristic 
of equilibrium using a sufficiently 'mild' nonequilibrium 
protocol. Many one-component systems indeed assem- 
ble in quasiequilibrium [10] [TTJ [T3HT5] if they are not 
deeply supercooled [16] [17] or plagued by slow particle 
dynamics [18]. However, a substantial literature sug- 
gests that multicomponent self-assembly is susceptible to 
kinetic factors even under conditions of weak driving [T9l - 
[26] . For instance, the Stranski- Totomanow assumption 
breaks down in a particular case of simulated binary 
colloid nucleation, where sluggish inter-species mixing 
prevents nuclei from establishing compositional equilib- 
rium [27] [28]. Also, binary crystals have been observed in 
simulation and experiment to grow out of compositional 
equilibrium, even under conditions 'mild' enough to pro- 
duce morphologically ordered structures [24-26 . To de- 
scribe such assembly theoretically, one must account for 
dynamical processes that drive a system away from low- 
lying paths on the free energy landscape [19 - 23 , 28U3T]. 
On a practical level, one can ask under what conditions 
can precisely-defined multicomponent structures be self- 
assembled, if evolution even near a phase boundary leads 
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to a structure not characteristic of the equilibrium one? 

Here we study this question within a prototypical ex- 
ample of multicomponent self-assembly. We apply simu- 
lation and quantitatively accurate analytic theory to the 
fluctuating growth of a model lattice-based fiber built 
from red and blue blocks. We show that when com- 
positional correlations of the equilibrium structure are 
not equal to those of the randomly-mixed material inci- 
dent on the fiber, the former cannot be generated at a 
finite rate of growth. Moreover, structures that assemble 
even close to the phase boundary can be very different 
to equilibrium ones. This absence of a quasiequilibrium 
regime occurs despite the fact that fibers near the phase 
boundary grow in a 'quasireversible' manner, displaying 
many unbinding events. However, by exploiting dynamic 
scaling connecting equilibrium and nonequilibrium pa- 
rameter manifolds, defined structures can be generated 
without error arbitrarily far from equilibrium. The fail- 
ure of a widely-made assumption for perhaps the simplest 
example of compositionally inhomogeneous self-assembly 
confirms the need for the development of dynamical the- 
ories [23] [26] [30] I32U35] in order to describe the self- 
assembly of multicomponent structures, and challenges 
the idea that the equilibrium structure is the natural tar- 
get for the self-assembly of components of many types. 

Model. We consider a 1-dimensional stochastic growth 
process in which a fiber is built from red and blue blocks 
(Fig. [TJa)). We add blocks to the right-hand end of the 
fiber with rate c (concentration). Added blocks are blue 
with probability Pbiue, and red otherwise. We allow the 
rightmost block to detach from the fiber with rate e _/3e % 
which depends on the nature of the rightmost bond of the 
fiber. Nearest-neighboring blocks of the same color inter- 
act with energy — 1{ = — e s , while the red-blue interaction 
is — ed (we set /3 = 1 throughout). We implemented this 
stochastic process in simulations by removing the ter- 
minal block with probability ^remove = V (l + ce^ ei ), 
where i = s or d as appropriate, and otherwise adding a 
block to the fiber end. This process is not designed to 



2 




FIG. 1: Fiber structures depend sensitively upon prepara- 
tion conditions, (a) Schematic of fiber energetics and growth 
rates, (b) Domain lengths £ of dynamically-generated fibers 
of given energy scales e s (simulations: solid lines; theory: 
dashed lines) are less than the corresponding equilibrium do- 
main lengths £o (horizontal dashed lines) for all points past 
the phase boundary c = Co (vertical dashed lines). This break- 
down of the quasiequilibrium assumption reflects the conflict 
between compositional correlations of unassembled material 
and the equilibrium structure, and occurs despite the fact 
that fibers close to the phase boundary grow in a 'quasire- 
versible' way. (c) Simulation data can be collapsed onto a 
master curve, summarizing the continuous variation of the 
scaled domain length g(£) = £(£ — 2)/(£ — 1) with concentra- 
tion c. 



model any particular physical system, but is instead a 
simple fluctuating protocol that captures the key kinetic 
constraint imposed upon a 3d structure in 3d space: ma- 
terial can be removed only from the interface between 
a structure and its environment [26]. By considering a 
lattice model in which complications of morphology and 
the possibility of internal block rearrangements are sup- 
pressed, we can explore directly how this constraint af- 
fects pattern generation. 

The 'equilibrium' we consider is the one often imag- 
ined in theories of self-assembly, corresponding to min- 
imization of the free energy of a structure of fixed 
size [lOj [TTJ [131 US • Such an equilibrium is achieved by 
our dynamical protocol in the limit of a large number of 
binding and unbinding events per site [50] : a key question 
is how the structure generated at a finite rate of growth 
compares to the equilibrium one. The energetics of a 
fiber of fixed size is that of the Id Ising model [36j 157] , 
with Hamiltonian T~L = — J^2 i SiSiJ r \ — h^2 i Si. Here 
the spin variable Si = ±1 describes a blue (Si = 1) or 



red (Si = — 1) block; the coupling J = (e s — e&)/2 is the 
penalty for domain wall (red-blue bond) creation; and 
the magnetic field h = — ln(l/pbiue — l)/2 describes the 
bias for blue blocks over red ones. Here we add red and 
blue blocks with equal likelihood, i.e. Pb\ue — 1/2, equiv- 
alent to h = 0. In this case the equilibrium structure 
of a fiber consists of equal proportions of red and blue 
blocks arranged into domains whose lengths t occur with 
probability p (£) = (f - l) _1 exp [£ln(l - ^q 1 )}. Here 
£o = 1 + exp(2/3 J) is the mean domain length in equilib- 
rium. 

Results. We carried out dynamic simulations for the 
choice ed = 1 and a range of values of e s > (for 
which £o > 2). For each set of energetic parame- 
ters we considered a range of concentrations c. When 
c < Co = 2/(e /3es + e^ ed ) the fiber does not grow. When 
c = Co the drift velocity of the fiber, averaged over dis- 
tances greater than a typical domain length, is exactly 
zero. The fiber therefore 'grows' only by diffusion of 
its rightmost end. When c > Co, the fiber grows with 
nonzero drift velocity. We therefore consider the concen- 
tration Co to define the 'phase boundary' between non- 
assembly and assembly. We stopped dynamic simulations 
when a fiber of length L = 2.5 x 10 4 blocks was gener- 
ated. We performed 10 4 simulations for each concentra- 
tion considered, except at the phase boundary, where the 
diffusive growth of a fiber was slow; there, we generated 
about 200 fibers for each set of conditions. For each en- 
semble of fibers we measured the probability distribution 
of domain lengths p(£) = n(£)/ X^=i where n(£) is 
the number of occurrences of domain length £ across all 
simulations at given thermodynamic conditions. These 
distributions were always exponential (Fig . Sl| ) with a 
mean £ depending on all three parameters c, e s and e^. 
This mean is shown in Fig. [ljb) for two choices of e s . 
At the phase boundary the equilibrium structure is gen- 
erated. For all points past the phase boundary, the dy- 
namic domain length £ is less than the equilibrium one 
£o, despite the fact that fiber growth close to the phase 
boundary is highly fluctional, exhibiting many unbind- 
ing events (Fig. [TJb), snapshot at right). Moreover, fiber 
structures are exquisitely sensitive to preparation condi- 
tions, and change continuously with supersaturation. 

Analytic theory reveals that this sensitivity is an in- 
evitable consequence of the different compositional statis- 
tics of the equilibrium structure and unassembled mate- 
rial. Consider defect variables rji = SiSi+i, where r\i — \ 
describes a bulk (same-color) bond and rji = — 1 is a de- 
fect (unlike-color) one. Let qb = [2rji — 1] be the likelihood 
that a given bond is a bulk one, where the average [•] is 
taken over many realizations of the dynamics. Enumera- 
tion of basic microscopic processes implies a drift velocity 
for bulk domains of 

«buik = c/2-^e-^, (1) 
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and a drift velocity for the fiber 

^fiber = C - 0e _/3es - (1 



(2) 



At the phase boundary, where ^buik = ^fiber = 0, 
Eqns. and ([2| give the equilibrium bulk fraction 
0o = 1/(1 + e^ ed-es )), and the equilibrium concentration 
Co = 2/(e^ ed + e^ €s ). £o an d Co are given by horizontal 
and vertical dotted lines on Fig. [ljb) . Fiber dynamics 
can be solved by requiring ^buik/^fiber = yielding 



c/2 - 0e 



-/3e s 



(1 



-/3e d ' 



(3) 



This equation is straightforwardly solved (Eq. Slh for 
the domain length £ = 1/(1 — 0), and we plot this so- 
lution as grey dashed lines in Fig. [TJb). The agreement 
with simulation is good, confirming the sensitivity of fiber 
structure to method of preparation. Moreover, the struc- 
ture of Eq. ([3| reveals the origin of this sensitivity. Its 
large-c limit returns the domain length characteristic of 
random mixing, ^ = 2, which is in general unequal to 
the equilibrium domain length £ - At the phase bound- 
ary, the balance of c-dependent and -independent terms 
in Eq. (|3| is such that bulk domains are generated at a 
rate characteristic of equilibrium. For finite supersatura- 
tion, however, the statistics of the fiber begins to reflect 
the statistics of random mixing: Eq. ([3| can be expanded 
in small deviations 5c = c — Co from the phase boundary 
to yield 



go (jp - 2)(|o - 1) 
co Co(Co-2) + 2 



6c. 



(4) 



Thus fiber structure is a continuous function of concen- 
tration, and the dynamic correlation length is less than 
the equilibrium one for any concentration c > Cq. Equiv- 
alently, structures generated at finite growth rate always 
sit above the minimum of the free energy landscape (Fig. 

Comparison of simulation and analytic theory thus re- 
veals that a fiber's dynamically- generated structure can 
be understood by considering only the relative net rates 
of bulk domain generation and fiber elongation. Al- 
though the time a given site has in which to relax is 
governed by the number of times the fiber end diffuses 
back and forth across it, with sites far from the fiber end 
being 'locked in' [26 , the structures generated by our 
fluctuating simulation protocol are explained solely by 
the competition between the compositional statistics of 
equilibrium and random mixing, with the departure from 
quasiequilibrium following from this competition enacted 
at a finite growth rate. The rate of this departure, given 
to lowest order by Eq. Q, has practical consequences 
for the design of fiber patterns. For target structures 
mimicking conventional binary crystals, having £o ~ 1, 
the conventional 'perturbative' scheme of growing some 




FIG. 2: Defined structures can be grown in equilibrium and 
far from it. Along each contour, defined by Eq. identical 
structures can be generated by 'reversible' dynamic protocols 
(at the phase boundary), or nearly irreversible ones (far from 
the phase boundary), illustrated by the simulation snapshots. 
Inset: 'distance' (£o — £)/£o ^° equilibrium of structures gen- 
erated along the £ = 15 contour. 



small distance 5c from the phase boundary leads to a 
structure with few errors, whose compositional correla- 
tions are numerically close to the target one (Fig. S3). 
We note that simulated 3d binary crystals can be grown 
close to compositional equilibrium [26 . However, when 
the desired target possesses long range compositional cor- 
relations, growth even close to the phase boundary results 
in many errors (Fig. S4). Similar predominance of kinet- 
ics is seen in weakly-driven growth of segregated binary 
assemblies [241125] . 

The 'perturbative' method of self-assembly thus be- 
comes increasingly unattractive as the compositional cor- 
relations of the target structure become increasingly com- 
plex. Here we propose an alternative 'nonperturbative' 
approach. Our model's dynamics, Eq. (J3|, can be cast in 
the scaling form 



c/c = g(€o)/g(€), 



(5) 



where g(£) = £(£— 2)/(£— 1) is a scaled correlation length. 
This expression reveals that the results of dynamic sim- 
ulations can be collapsed onto a straight line connecting 
a manifold of nonequilibrium conditions with an equilib- 
rium point at (1,1). We verified this collapse for a range 
of simulations carried out at 40 different phase points 
(Fig. []Jc)) whose values of e s ranged from 2 to 6. The 
nature of the master curve demonstrates that no fiber 
having £ = £o can be generated anywhere past the phase 
boundary c = Co, and emphasizes the sensitivity of struc- 
ture £ to preparation condition c. But it also describes 
a useful connection between patterns in and out of equi- 
librium. It can be rearranged to read 



_ 2(g- 1) , o €d _ * e v 



(6) 



revealing that structures of identical £ can be generated 
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on a manifold of parameters c, e s and ej, regardless of dis- 
tance to the phase boundary. We illustrate this manifold 
in Fig. |2j £ is constant along the displayed contour lines. 
Simulations verify this prediction (Fig. S5): structures 
generated dynamically at the circles on each contour line 
are indistinguishable. However, fiber generation changes 
from being purely diffusive (and very slow) at the phase 
boundary, to being nearly irreversible (and very fast) far 
from it, illustrated by the space-time trajectories shown 
(Fig. S5). Moving rightwards along contours, the gener- 
ated structures lie increasingly far from the equilibrium 
one for the corresponding energetic parameters (Fig. [2] 
inset); nonetheless, structures characteristic of equilib- 
rium at a particular point on the phase boundary can be 
generated by a continuum of nonequilibrium protocols 
away from it. In this sense, equilibrium structures can 
be grown, error-free, arbitrarily far from equilibrium. 

Conclusions. Compositionally inhomogeneous struc- 
tures are found abundantly in biology [38j [39] and are 
increasingly the target of designed self-assembly [25j l40l 
[44] . We have shown that the growth of a two-component 
fiber whose compositional correlations are not random 
never results in the equilibrium structure, recapitulating 
the importance of kinetic factors to the self-assembly of 
multicomponent structures identified in previous theo- 
retical [T9H23] . simulation [24H27] and experimental [25] 
work. It is increasingly clear that for multicompo- 
nent self-assembly, specification of a free energy sur- 
face [45j [46] is in general not enough to ensure assem- 
bly of the desired structure, even under conditions of 
weak driving. Instead, explicit accounting of how micro- 
scopic dynamics [47 select assembly pathways [481 EH] 
is required. Indeed, we have shown that by exploiting 
dynamic scaling of the rates of basic microscopic pro- 
cesses, precisely-defined structures can be self-assembled 
out of equilibrium. Two questions arise naturally from 
this work: 1) Which classes of real multicomponent ma- 
terials can be grown in quasiequilibrum? And 2) for those 
that cannot, is the 'nonperturbative' self-assembly strat- 
egy described here viable? We speculate that the answer 
to question 1) is 'those structures displaying local free 
energy barriers on the scale of the desired compositional 
correlation length'. 3d binary crystals with different lat- 
tice structures appear to display different propensities for 
growing in compositional equilibrium [26] . and we spec- 
ulate that the free energy cost of adding a few particles 
to an existing structure may be a key difference in this 
regard. However, when the desired compositional corre- 
lation length is large, achieving free energy barriers over 
the extent of this length is likely unworkable. We then 
turn to question 2), which experimental evidence hints 
might be answered in the positive. The dynamic scaling 
that permits the far-from-equilibrium assembly strategy 
we have described is strikingly reminiscent of data col- 
lapse seen in segregated binary structure growth, where 
an effective structural order parameter scales with effec- 



tive crystal growth velocity [25]. Such scaling, occurring 
in a regime in which crystals are morphologically ordered, 
suggests the possibility of designing compositionally or- 
dered structures in equilibrium, and assembling them far 
from it. 
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